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The leading-order parton processes that produce a dilepton with large transverse momentum 
predict that the transverse polarization should increase with the transverse momentum for almost 
any choice of the quantization axis for the spin of the virtual photon. The rate of approach to 
complete transverse polarization depends on the choice of spin quantization axis. We propose axes 
that optimize that rate of approach. They are determined by the momentum of the dilepton and 
the direction of the jet that provides most of the balancing transverse momentum. 
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The spins of particles produced in high energy colli- 
sions carry important information about the fundamen- 
tal interactions between elementary particles, but it is 
difficult to access that information. Since the spin can- 
not be measured directly, information about the spin of 
a particle must be inferred from the angular distribution 
of its decay products. The accessibility of information 
about the spin is affected by the choice of a spin quan- 
tization axis (SQA), because some angular variables are 
more easily measured than others. Hadron collisions pro- 
vide an additional complication, because the fundamen- 
tal interactions involve collisions of partons with varying 
longitudinal momenta. Integration over these longitu- 
dinal momenta tends to dilute the information carried 
by the spins of the final-state particles. This raises an 
important question: which SQA will maximize that in- 
formation? 

The simplest process for which this question can be ad- 
dressed is the production of a dilepton, a lepton and an- 
tilepton with opposite charges. There are various sources 
of dileptons in high energy collisions, but we focus on pro- 
duction through a virtual photon, which can be regarded 
as a spin-1 particle with a variable mass equal to the 
invariant mass of the lepton pair. There are some sim- 
ple dilepton production mechanisms that predict that the 
virtual photon should be transversely polarized. For the 
Drell- Yan mechanism, the annihilation of a quark and an 
antiquark into a virtual photon, the vector coupling to 
the quark and antiquark ensures that the virtual photon 
is transversely polarized 0]. The leading-order parton 
processes for a dilepton with large transverse momentum 
are dominated asymptotically by photon fragmentation, 
the production of a transversely-polarized real photon 
followed by its decay into a collinear lepton pair [2j. For 
these examples, the leading-order prediction of pertur- 
bative QCD is that the virtual photon is transversely 
polarized. Longitudinal virtual photons are produced 
by higher-order processes at a rate that depends on the 
SQA. A reasonable criterion for an optimal SQA is that it 
minimizes the cross section for producing a longitudinal 
virtual photon. 



For the Drell- Yan mechanism, an optimal SQA was 
identified long ago: it is the Collins-Soper axis 0]. In 
this paper, we derive optimal SQA's for the parton pro- 
cesses that create a dilepton with large transverse mo- 
mentum. We apply them to dilepton production at Fer- 
milab's Tevatron and at CERN's Large Hadron Collider 
(LHC). 

We consider the production of a dilepton whose invari- 
ant mass Q is well below the Z° resonance. The orienta- 
tion of the lepton momentum in the dilepton rest frame is 
given by a polar angle 8 with respect to a chosen axis and 
an azimuthal angle <f> with respect to a chosen plane con- 
taining that axis. If <fi is integrated over, the differential 
cross section can be expressed as 

da a fdQ 2 \ I + x 2 ,„ 2 J 

where x = cos 9 and o~t and ctl are the cross sections for 
a transverse and a longitudinal virtual photon, respec- 
tively. The polar axis that defines the angle 9 can be 
identified with the SQA of the virtual photon. Our ques- 
tion can then be restated: for which SQA is the virtual 
photon most strongly polarized? 

The longitudinal polarization 4- vector er, for a virtual 
photon with 4-momentum Q must satisfy Q ■ el = and 
e\ = — 1. The most general 4- vector satisfying these 
conditions can be written in the form 

e J=l7V^, X» = {-g^ + Q»Q V IQ 2 )X V , (2) 

where A is a 4-vector. The physical interpretation of X 
is that in the dilepton rest frame, which is defined by 
Q = 0, the direction of the 3- vector —X is the SQA of 
the virtual photon. 

We now consider the production of a dilepton in the 
collision of two hadrons with 4-momenta Pi and P 2 . The 
SQA is generally chosen to lie in the production plane 
defined by the momenta of the colliding hadrons in the 
dilepton rest frame. Thus X in Eq. (J2]) has the form 

X" = aPf + &P 2 ", (3) 
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where a and b are scalar functions. In the dilepton rest 
frame, the SQA is antiparallel to the unit vector X = 
(a Pi + bP 2 )/\aPi + bP 2 \, so it is determined by the ratio 
a/b. If the leptons are the only particles in the final state 
whose momenta are measured, then a/b can only depend 
on the 4- vectors Pi, P 2 , and Q. If additional information 
about the final state is measured, a/b can also depend on 
this information. 

If the invariant mass Q = (Q 2 ) 1 ^ 2 of the dilepton is 
large compared to the scale Aqcd of nonperturbative 
effects in QCD, the inclusive cross section can be cal- 
culated using QCD factorization formulas. Suppose the 
cross section is dominated by a specific parton process 
which, at leading order in the QCD coupling constant, 
predicts that the virtual photon is transversely polarized 
for any choice of the SQA. Then a natural prescription 
for an optimal SQA is that it minimizes the longitudi- 
nal cross section from that parton subprocess at leading 
order. Such an optimal SQA will depend on the longitu- 
dinal momentum fractions of the colliding partons, which 
cannot be directly measured. 

We first consider the production of a dilepton with 
large invariant mass Q 3> Aqcd- The Drell-Yan mecha- 
nism is the parton process qq — ► 7*. If the momenta of 
the quark and antiquark are collinear to those of their 
parent hadrons, the virtual photon is transversely polar- 
ized for any choice of the SQA. A rigorous QCD calcu- 
lation of the transverse momentum distribution for the 
dilepton requires the resummation of the effects of the 
emission of soft gluons from the colliding partons Q . One 
of the most important qualitative effects of soft-gluon 
emission is that it gives transverse momenta to the collid- 
ing partons. A simple model for these effects is the parton 
model with intrinsic transverse momentum, in which the 
momentum distribution of a parton is the product of the 
parton distribution and a function of the 2-dimensional 
transverse momentum vector . Collins and Soper cal- 
culated the angular distribution of the dilepton in this 
model [3|. The cross section for qq — » 7|_, where 7^ is a 
longitudinal virtual photon, is 

. 8n 2 e 2 q a(k 2 ± )(a 2 x 2 + b 2 x 2 ) sf 

aL= ZQ 2 (ax 2 -bxi) 2 6 ^ S ~Q)> ( 4 ) 

where xi and x 2 are the longitudinal momentum frac- 
tions of the colliding partons and e q is the electric charge 
of the quark. We have used the expression for the lon- 
gitudinal polarization vector in Eqs. ((2]) and (|3|). The 
cross section has been expanded to second order in fcu_ 
and k 2 ± and averaged over them using (k l nJ _) = and 
(k l n j_k 3 n j_) = \{k 2 _)5 % /[, where S l f is the unit tensor in the 
transverse dimensions. 

Collins and Soper proposed a convenient set of angles 
9 and d> for the lepton momentum in the dilepton rest 
frame [3|. Lam and Tung pointed out that the Collins- 
Soper axis is an optimal SQA for the Drell-Yan mecha- 
nism Q. It minimizes the cross section for longitudinal 



virtual photons in Eq. ((4|) , thus making the leading-order 
prediction that the virtual photon will be transversely 
polarized as robust as possible with respect to radiative 
corrections that generate transverse momenta for the col- 
liding partons. Minimizing the cross section in Eq. ^ 
with respect to the ratio a/b, we find 



a/b\ 



as 



-xi/x 2 . 



(5) 



Under the assumption that the cross section is dominated 
by the Drell-Yan mechanism qq — > 7*, we can derive an 
expression for xi/x 2 in terms of variables that can be 
directly measured. At leading order in fc„j_, the energy- 
momentum conservation condition is Q — xiPi + x 2 P 2 , 
which implies 



x 1 /x 2 \ 



cs 



Q-P 2 /QPi. 



(6) 



Equations (O and (J6j) define the Collins-Soper axis. The 
corresponding 4-vector X in Eq. becomes 



A cs - 



pfl 



Q-Pi 



2P 2 



(7) 



In the dilepton rest frame, the SQA is along Vi — v 2 , 
where the v^s are the velocities of the colliding hadrons. 

We now consider the production of a dilepton with 
large transverse momentum Qt ^> Aqcd ■ At leading or- 
der in a s , the parton processes that create a virtual pho- 
ton at large Qt are qq — > 7*5, qg —> 7*5, and qg — > 7*9. 
For Qt ^> Q, the cross section for a longitudinal vir- 
tual photon is suppressed by a factor of Q 2 /Q T . The 
coefficient of Q 2 /Q T depends on the choice of SQA. Our 
prescription for the optimal SQA is to minimize the lon- 
gitudinal cross section, which makes the transverse po- 
larization for fixed Q 2 increase as rapidly as possible with 
increasing Qt- We will refer to the optimal SQA for the 
parton process ij — > 7*fc as the optimal ij axis. 

We first determine the optimal qq axis. At leading or- 
der in a s , the dependence of the differential cross section 
for qq — * ^* L g on a and b is 



a°Q (awi 



a 2 x 2 + b 2 x\ 



bwo 



abQ 2 s ' 



(8) 



where wi = Q ■ Pi, w% = Q • P2, and s is the center-of- 
momentum (cm.) energy of the colliding hadrons. There 
is a delta function constraint on these variables: 2(xiwi + 
x 2 w 2 ) = x\x 2 s + Q 2 . Minimizing Eq. (jHJ) with respect to 
a/b, we get 



xlwl + Z 1 ' 2 



x\{2wiw 2 - Q 2 s) 



(9) 



Z = {x\w\ + x\w\) 2 — x\x\Q 2 ''s(AwiW2 — Q 2 s). 



This also gives the optimal qq axis. 

We next determine the optimal qg axis. At leading or- 
der in a s , the dependence of the differential cross section 



3 



for qg — > ^* L q on a and b is 

(ax 2 - bxi) 2 + b 2 x\ 



(PQ (aw\ + bvj2) 2 — abQ 2 s 
Minimizing with respect to a/b, we get 

2x\w\ - x\w\ + Z 1 ' 2 



(10) 



(11) 



^7* 9 X 2 (2xiwl + 2x 2 WiW 2 - x 2 Q 2 s 

Z = (2x 2 w 2 + x 2 w 2 + 2xix 2 wiiu 2 — x\x 2 Q 2 s) 2 
-x\x\Q 2 s(4:W 1 w 2 - Q 2 s). 

This also gives the optimal qg axis. The values of a/b for 
the optimal gq and gq axes are obtained by taking the 
reciprocal of the expression on the right side of Eq. (TTTj) 
and then interchanging x± and w\ with x 2 and W2- 

The expressions for a/b in Eqs. © and (TTTj) depend 
on the longitudinal momentum fractions x\ and x 2 of 
the colliding partons only through the ratio x\/x 2 . Our 
optimality criterion was based on the assumption that 
a specific parton process dominates. If that parton pro- 
cess implies a value for X\ /x 2 that depends on a measur- 
able property of the final state, we can insert that value 
into Eqs. ([9]) and (fTTj) to obtain optimal SQA's that are 
experimentally useful. In the leading-order parton pro- 
cesses for a dilepton with large Qt, the large transverse 
momentum is balanced by that of a single parton. This 
recoiling parton produces a jet of hadrons whose mo- 
menta are nearly collinear to that of the parton. The 
polar angle of the recoiling parton in the hadron cm. 
frame is approximately equal to the polar angle 0j e t of 
the jet. The ratio X\/x 2 can be expressed as a function 
of #jet and the transverse and longitudinal momenta Qt 
and Ql of the dilepton in the hadron cm. frame: 



■i'2 



(Ry* + Q L ) sinfljet + QtQ- + cosfr 



jet J 



optimal (i? 7 * — Ql ) sin I 



jet 



Q T (l - cos6> jct ) ' 

(12) 

where E 7 * = (Q T + Q\ + Q 2 ) 1 / 2 . Inserting this into 



Eqs. j9| and (fTT|) . we obtain expressions for a/b for the 
optimal qq and qg axes that depend on quantities that 
can be directly measured. 

Beyond leading order in a s , there can be more than 
one parton in the final state, so there can be more than 
one jet with large transverse momentum. As a general 
prescription for x\/x 2 , we choose 6*j c t in Eq. (|T2]) to be 
the angle in the hadron cm. frame of the jet with the 
largest transverse energy. Since the direction of the jet is 
insensitive to soft gluon radiation and to the splitting of a 
parton into collinear partons, QCD radiative corrections 
to the polar angle distributions defined by our optimal 
SQA's can be calculated systematically using perturba- 
tive QCD. The number of jets in an event depends on 
the size of the angular cone used to define the jet. If a 
reasonably large value is chosen for the cone size, most of 
the events will contain a single jet with large transverse 



momentum balancing that of the dilepton. A small frac- 
tion of the events will have more than one such jet. For 
these multijet events, the expression on the right side of 
Eq. (|12p may not be a good estimate for x\/x 2 . However 
since these events are only a small fraction of the total 
number of events, they should not seriously decrease the 
transverse polarization. Our expression for x\/x 2 is not 
useful for fixed-target experiments, because the recoiling 
jet is usually not observed. However #j 0t can be measured 
relatively easily in high energy hadron colliders. 

To illustrate the use of our optimal SQA's, we ap- 
ply them to the production of a dilepton with invari- 
ant mass Q = 3 GeV at the Tevatron, which is a pp 
collider with cm. energy 1.96 TeV. We calculate the 
cross sections for a longitudinal virtual photon using the 
CTEQ6L parton distributions Q with 3 flavors (u, d, and 
s). The factorization and renormalization scales are set 
to fi = (Q 2 + Q\) x l 2 ■ For the strong coupling constant 
a s (fi), we use the next-to-leading order formula with 4 
flavors of quarks and Aqqd=326 MeV. For the electro- 
magnetic coupling constant a, we use 1/132.6. We im- 
pose a rapidity cut \y\ < 1 on the dilepton momentum. 

We compare our optimal SQA's with three other 
choices for the SQA: 

(i) the cm. helicity axis defined by X cm h = Pi + P 2 . 
The projection of the spin of the virtual photon 
along this axis is minus its helicity in the hadron 
cm. frame. 

(ii) the Collins-Soper axis defined by Eq. ([7]). 

(iii) the perpendicular helicity axis defined by 



r 2 



Q-Pi Q-P 2 



(13) 



In the dilepton rest frame, the perpendicuclar helicity 
axis is along Vi + v 2 . The spin of the virtual photon 
along this axis coincides with its helicity in the frame 
obtained from the hadron cm. by a longitudinal boost 
that makes the dilepton momentum perpendicular to the 
beam direction. The perpendicular helicity axis coincides 
with the cm. helicity axis when the rapidity of the dilep- 
ton is 0. One advantage shared by X±h and Acs is that 
they are invariant under independent longitudinal boosts 
of the two colliding hadrons. Thus they can be expressed 
equally well in terms of the colliding parton momenta or 
the colliding hadron momenta. 

A convenient polarization variable for a virtual photon 
is a = (ot — 2ol)/(pt + 2<jl)- It can be measured from 
the polar angle distribution of the lepton, which is pro- 
portional to 1 -I- a cos 2 8. The leading-order predictions 
for a for various choices of SQA are shown as functions of 
Qt/Q in Fig. [TJ For the Collins-Soper axis, a is negative 
for Qt > Q- For the cm. helicity axis, the perpendicular 
helicity axis, and our two optimal axes, a approaches 1 
at large Qt, indicating pure transverse polarization. The 
approach to 1 is significantly faster for the perpendicular 
helicity axis than for the cm. helicity axis. However it 
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FIG. 1: (color online). Polarization variable a for various 
spin-quantization axes as functions of Qt/Q for dileptons 
with invariant mass Q — 3 GeV at the Tevatron. 

is even faster for both of the optimal axes. The two op- 
timal axes give essentially the same a for Qt > 2Q, but 
for Qt < 2Q a is closer to 1 for the optimal qg axis than 
for the optimal qq axis. 

The requirement that the recoiling jet be captured by 
the hadronic calorimeter imposes a constraint on its pseu- 
dorapidity rjj e t = lntan(0j C t/2). To provide some idea of 
how much this additional constraint might decrease the 
data sample, we impose the pseudorapidity cut |^j e t| < 1- 
For dileptons with Q — 3 GeV at the Tevatron, the frac- 
tion of events satisfying the dilepton rapidity cut that 
also survive the jet pseudorapidity cut is greater than 
0.5 for Qt > Q- This fraction is large enough that mea- 
suring the angle (9j Ct should not dramatically decrease the 
size of the data sample. 

We have also compared the various quantization axes 
for dilepton production at the LHC, which is a pp col- 
lider with cm. energy 14 TeV. We imposed a rapidity cut 
\y\ < 3 on the dilepton. For Qt > Q, the leading-order 
results for a differ by less than 0.01 from the results for 
the Tevatron in Fig. [T] for all the SQA's with one excep- 
tion. For the cm. helicity axis, a at the LHC is smaller 
by more than 0.1 in the region near Qt — 1.5 Q. The 
leading-order predictions for a as functions of Qt/Q are 
similar to those for the Tevatron in Fig. [T] If we impose 
a pseudorapidity cut |?7j c t| < 3 on the recoiling jet, the 
fraction of events satisfying the dilepton rapidity cut that 
also survive the jet pseudorapidity cut is greater than 0.8 
for Q T > Q. 

The CDF and DO Collaborations have measured a as 
a function of Qt for dimuons produced at the Tevatron 
from decays of charmonium mesons 0] and bottomonium 
mesons [8j, |9j. The variable a for the cm. helicity axis 



was measured for Qt/Q as high as 9.7 for J/ip and 2.1 for 
T(15). The data samples collected by CDF and DO are 
large enough that it should also be possible to measure a 
for the sidebands of these quarkonium resonances, which 
come from dileptons produced by virtual photons. At 
the LHC, it should be possible to measure a out to much 
larger values of Qt/Q- 

Our leading-order results for a as a function of Qt in- 
dicate that the virtual photon will be significantly more 
strongly polarized at large transverse momentum if we 
use optimal SQA's. A quantitative prediction of the po- 
larization requires calculations to next-to-leading order 
(NLO) in a s . The angular distributions for lepton pairs 
at large Qt have been calculated to NLO by Mirkes and 
Ohnemus [lj| • At NLO, there can be two partons in the 
final state. In Ref. [T3|, the momenta of the two partons 
were integrated over, so they cannot be resolved into two 
separate jets. In order to use these NLO results, we would 
have to modify our prescription for x\jxi- At small Qt, 
it is also necessary to sum up the effects of soft-gluon 
emission to all orders 

Similar methods could be used to derive optimal SQA's 
for the production at large transverse momentum of mas- 
sive particles, such as J/ij), T, and the weak bosons W ± 
and Z°. With an optimal SQA, these states should be 
more strongly polarized, so measurements of a will pro- 
vide more information. It may also make the theoretical 
predictions more robust with respect to radiative correc- 
tions. Optimal SQA's might be useful for determining 
the spin of new particles created at the LHC. They pro- 
vide a new window into the spins of particles created by 
parton collisions. 
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